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Abstract 

We study the addition of flavor to the gravity dual of = 4 three-dimensional gauge 
theories obtained by wrapping D4-branes on a two-cycle of a non-compact Calabi-Yau 
two-fold. In this setup the flavor is introduced by adding another set of D4-branes that 
are extended along the non-compact directions of the Calabi-Yau which are normal to the 
cycle which the color branes wrap. The analysis is performed both in the quenched and 
unquenched approximations. In this latter case we compute the backreacted metric and 
we show that it reproduces the running of the gauge coupling. The meson spectrum and 
the behavior of Wilson loops are also discussed and the holographic realization of the Higgs 
branch is analyzed. Other aspects of this system studied are the entanglement entropy and 
the non-relativistic version of our backgrounds. 
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1 Introduction 



The gauge/gravity correspondence [T] is one of the major achievements of string theory in 
the last ten years (see |2j for a review). Its extension to include more realistic theories is, 
clearly, a topic of great interest. In particular, the addition of matter degrees of freedom is 
essential to get nearer to the string theory description of QCD-like theories. This addition 
can be performed by including extra (flavor) branes [3], which should be extended along all 
the gauge theory directions and wrap a non-compact cycle in order to make its worldvolume 
symmetry a global symmetry from the point of view of the gauge theory living in the color 
branes. They introduce an open string sector that corresponds to having hypermultiplets 
transforming in the fundamental representation in the gauge theory side. 

If the number of flavors, Nf, is small compared with the number of colors, Nc, one can 
treat the flavor branes as probes in the supergravity background created by the color branes. 
This deflnes the so-called quenched approximation, which corresponds, in the fleld theory 
side, to neglecting quark loops, which are suppressed by factors 1/Ncin the 't Hooft large Nc 
expansion [1]. This quenched holographic flavor has been explored extensively in the past 
few years. In particular, by analyzing the normalizable fluctuations of the probe branes, the 
spectra of mesonic excitations of different theories have been analyzed (see [5] for a review 
and a list of references). 

On the contrary, if the number of flavors is of the order of the number of colors {Nf ~ Nc) 
the backreaction of the flavor branes on the geometry cannot be ignored and one has to deal 
with a system of gravity plus branes, the latter acting as dynamical sources for the different 
supergravity flelds. On the fleld theory side the inclusion of the backreaction in this Nf ~ Nc 
regime is equivalent to considering the so-called Veneziano limit [6] , in which Nc and Nf are 
large and their ratio Nf/Nc is flxed. In this limit quark loops are not suppressed and the 
flavor is unquenched. 

Another direction in which the original gauge/gravity duality has been generalized is by 
extending it to theories with lower amounts of supersymmetry. A general strategy to carry 
out this extension of the correspondence to less supersymmetric models is to consider higher 
dimensional branes wrapping cycles. At energies small compared with the size of the cycle 
the theory becomes effectively four-dimensional. Moreover, the gauge theory living on the 
worldvolume of the wrapped brane has to be topologically twisted in order to preserve some 
fraction of supersymmetry. Examples of duals of Ad gauge theories constructed in this way 
are the geometries analyzed in for Af = 1 and in [8] for A/" = 2, which correspond to 
D5-branes wrapping a two-cycle inside a Calabi-Yau (CY) manifold. 

In this paper we will study the dual of A/" = 4 gauge theories in three space-time dimen- 
sions, obtained by wrapping D4-branes of the type IIA theory on a two-cycle of a Calabi-Yau 
two-fold. The corresponding unflavored supergravity solution was found in [9j and studied 
in detail in [10], where it was shown to reproduce the exact perturbative running coupling 
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constant and the metric of the moduh space of the gauge theory. In the present paper we will 
analyze the addition of flavor branes in this setup. These flavor branes are also D4-branes, 
which are extended along non-compact directions of the CY two-fold in such a way that no 
further supersymmetry is broken. 

We will analyze the addition of flavor to the 3d, Af = A theory both in the quenched 
and unquenched approaches. In the quenched formalism we will compute numerically the 
spectrum of mesonic excitations and we will be able to find a very good analytical estimate 
of the masses. We will also study the Wilson loops and the corresponding quark- antiquark 
potentials. In this latter case we will discover that a phase transition is produced when some 
of the parameters of the solution are varied and we will compute the corresponding critical 
exponents. 

We will include the backreaction by following the proposal of ref. [H], in which the 
localized brane sources are substituted by a continuous distribution of flavor branes in their 
transverse directions (see also for a similar analysis in the context of non-critical string 
theory). This smearing brane approach has been applied successfully to various string duals 
in several dimensions ([l3]-[25]). It has the advantage that it avoids having the 5- function 
sources of the localized branes in the BPS equations and, therefore, makes searching for 
solutions more feasible. 

In our case we will start with the same ansatz for the metric as in the unflavored system. 
However, the presence of flavor branes modifies the Bianchi identity of the RR field strength 
F4, which forces us to modify the unflavored ansatz for F4. Once this fact is taken into 
account, one can get a system of first-order differential equations by imposing that the 
supersymmetric variation of the gravitino and dilatino of type IIA supergravity vanish. These 
first-order equations are simple, but they are difficult to solve in general. However, we will 
see that they can be solved analytically in certain regions of the space. Remarkably, this 
is enough to compute the modification of the running of the coupling constant due to the 
matter hypermultiplets in the gravity solution. We will show that this modification matches 
exactly the field theory results. 

The organization of the rest of this paper is the following. We will start in section [2] 
by reviewing the unfiavored background. We will set up our notations and write the BPS 
equations and their solution for the unfiavored system. In section [3] we study the addition of 
fiavor branes. In particular, we obtain the BPS equations for the backreacted background, 
which are then numerically integrated. Section H] is devoted to analyzing the matching 
between our fiavored background and its field theory dual. By means of a probe calculation 
we verify that the running of the gauge coupling with fiavor is reproduced by our solution. 
We also discuss in this section how to realize the Higgs branch in our setup. 

The analysis of the meson spectrum is the object of section [51 We study the mass levels 
both in the quenched and unquenched solutions and we also consider the meson spectra 
in the Higgs branch. In section [6] we explore the behavior of the Wilson loops. In our 
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study of the energy of the quark-antiquark pair we will find some critical phenomena and 
we will evaluate the corresponding critical exponents. Section [7] is devoted to presenting our 
conclusions and summarizing our results. 

The paper is completed with several appendices. In appendix |X] we give details of the 
derivation of the BPS equations for the general flavored system and we check that the 
equations of motion of the system are satisfied if the BPS equations hold. In appendix [Bl we 
find the supersymmetric embeddings of D4-brane probes in our background. In appendix O 
we find additional solutions of the unflavored equations and we obtain the background dual 
to a non-relativistic system that can be generated by deforming our solutions. Finally, in 
appendix O we analyze the entanglement entropy for our model in the UV. 



2 The unflavored background 



Following the analysis of refs. [9j and [lO], let us consider the background of type IIA 
supergravity created by a stack of Nc D4-branes wrapped on a two-cycle C2 of a Calabi-Yau 
cone of complex dimension two, according to the following brane setup 
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where S"^ represents the directions of the two-cycle (which is topologically a two-sphere) and 
N2 are the directions of the normal bundle to C2 ■ In the above setup a circle represents 
wrapped directions, whereas the symbols "-" and "." denote unwrapped worldvolume and 
transverse directions respectively. We shall parameterize the cycle C2 by means of two 
angular coordinates {9, 0) with < ^ < vr and < < 27r and we will denote by a the radial 
coordinate of the CY cone which, together with an angular coordinate ip, will parameterize 
the normal bundle A^2- Moreover, we shall choose a system of spherical coordinates for the 
transverse M^, r being the corresponding radial coordinate and {6, (p) the angular variables 
(0<^^<7r, O<0< 27r). The concrete ansatz for the ten-dimensional string frame metric 
we will adopt is the following: 



2<S> 



+ e 



-2$ 



dx'i ^ + [ dO' + sin^ M^'^ 

_ _\ 2 



^ f c/a^ + (T^ [dip + cos dd^y\ + dr'^ + [dd'^ + sin^ Odcp^) 



{21] 



where dx\ 2 denotes the Minkowski metric in 2+1 dimensions and the range of is < < 
27r. Notice that %p is fibered over the {9, (p) two-sphere. For convenience we have included in 
(12. ip the radius i?, given by: 

= STxg,N,{a')^ , (2.2) 
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with gs and a' being respectively the string couphng constant and the Regge slope. The 
ansatz (12.11) contains two functions: Z, which controls the size of the cycle, and which is 
the dilaton of the type IIA theory. Both of them should be considered as functions of the 
two radial coordinates r and cr: 

$ = $(r,a) , Z = Z{r,a) . (2.3) 

As in any other background generated by D4-branes, the ansatz should be endowed with an 
RR four-form F4. Let C3 denote the corresponding three-form potential (F4 = dC^). We 
shall adopt the following ansatz for C3: 

C3 = - guj2 A (di) + cos 9 d(j)) , (2.4) 

where g{r, a) is a new function and uj2 is the volume element of the {6, 0) two-sphere: 

UJ2 = sinOde A d(j) . (2.5) 

The corresponding RR four-form field strength will be: 

-F4 = — {gda + g dr) /\ UI2 /\ {dip + cos9d(p) + g 0J2 A U2 , (2.6) 

where we have denoted: 

' = dr, ■ = d^. (2.7) 

We will require that our background preserves eight supersymmetries. which is the appro- 
priate number of SUSYs for a supergravity dual of A/" = 4 gauge theories in three dimensions. 
As shown in detail in appendix lAj the vanishing of the different components of the super- 
symmetric variations of the gravitino and dilatino gives rise to a system of first-order BPS 
equations for the functions $, Z and g entering our ansatz. This system is the following: 

g = -R^ Z' , e-^* a = ZZ , 

g' = _4ar2e-^*$ , g = -a R~^ Z-^ e"^"^ g + 4 a e"^* $' . (2.8) 

It is interesting to notice that not all the equations in (12. 8p are independent. Actually, one 
can check that the equation for g in (12.81) can be obtained from the others. Moreover, one can 
combine the different equations in (12.81) and get a single second-order PDE for the function 
Z{r, a), namely: 

rZ(^Z-aZ^ = a i^r Z'^ + r Z" + 2Z'^ . (2.9) 

Notice that, if Z is known, the other functions $ and g can be determined from the first 
two equations in (12.81) . Moreover, we check in appendix [A] that the second order equations 
of motion for the RR four- form F4, dilaton $ and the metric Gmn of type IIA supergravity 
follow from the system (12.81) . 
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2.1 Integration of the BPS system 



The BPS system (12.81) can be integrated by elementary methods when a = and r varies. 
Indeed, it follows from the last line in (12.81) that g{r, cr = 0) is constant. Let us put: 

g{r,a = 0)=go . (2.10) 
Then, the first equation in (12.81) for cr = can be readily integrated, namely: 

Z'{r,0) = - ^ ^ Zir,0) = + constant. (2.11) 

The actual value of Qq can be obtained from the quantization condition of the RR four-form 
flux F4. Actually, the best way to perform this analysis is by using the approach in which 
the solution is obtained by uplifting from gauged supergravity. This study was done in 
refs. [9[ [To] and allows one to find a solution of the system (12.81) for arbitrary values of the 
variables r and a. Here we will just reproduce this solution with our notations. First of all, 
let us define the function T{Z) as follows: 

1 OK" 

T{Z) =Z, + 2{Z -Z,) + —{Z- Z,f , (2.12) 

Zjf 

where Z^, and k are constants. Then, Z{r, a) is determined implicitly as the solution of the 
equation: 



[Z^ - Z) 



2 (y^ R 



+ -^-^ " = - . (2.13) 
T{Z) J 8 ^ ' 

We are able to solve this quartic equation in Z exactly, although we must be careful that 
when we use this solution in the following we are always picking up the appropriate root. For 
each value of r and cr it must be checked that the solution is real and less than Z^,. Notice 
from (12.131) that Z^ is just the constant value approached by Z in the UV, i.e. when either 
r or cr is large. Moreover, g{r, a) and $(r, a) are given by: 



^ R^Z Z,)r^ 4^ ^ 2 Z{Z, Z 

^ r'i A —(7"^ r „ 1 \ ) 



r2 J ^rr2 



T{Z) 



As a check of this solution, one can easily prove by taking derivatives of (I2.12l) - (l2.13p that 
the BPS system (12. 8p is satisfied. Moreover, by taking cr = in (I2.13P one can verify that 
Z{r, 0) is indeed of the form (12. lip , with go being given by: 

9^=-^- (2-15) 
Using this result we can rewrite Z{r, 0) as: 

Z(r,0) = - ^ . (2.16) 
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Furthermore, by using (12.161) to evaluate the right-hand side of the second equation in fl2.14p . 
one can obtain the value of the dilaton at a = 0, namely: 



7 — K r3 
^-4*M) ^ tL . (2.17) 



From the explicit expressions for Z and $ written above one easily concludes that, when 
r and a are small enough, the supergravity solution is not valid because the function Z 
becomes negative and/or the dilaton $ becomes complex. This phenomenon is related to 
the so-called enhangon mechanism [26] (see below). We can estimate the scale at which this 
mechanism occurs by computing from (I2.16P the value of r for which 2(r, 0) vanishes. This 
determines the so-called enhangon radius re, given by: 

(2.18) 



8Z 



Notice also that the sign of the right-hand side of (I2.17P becomes negative for sufficiently 
small r. Actually, the numerator in (I2.17P changes its sign precisely at r = rg , whereas 
the change of sign of the denominator depends on the value of the constant k. Indeed, by 
analyzing the discriminant of the quadratic function in the denominator of (12.171) one easily 
concludes that this equation has no real roots for k > 1/16. Thus in this case e"^**-'"''^'' 
changes its sign precisely at the enhangon radius Vg. However, for k < ^ this change of sign 
occurs for larger values of r, namely for r = Vh, where is given by: 



rn = {I + Vl-16/t) re . (2.19) 
Clearly, when k < 1/16 the space ends, at a = 0, when r = th > rg. 

2.2 UV form of the metric 

As mentioned above, it follows from (12.131) that the function Z approaches the constant value 
Z = Z^ a,s one moves into the UV region. Since the solution of the algebraic equation (12.130 
for 2(r, cr) is complicated, we can try to perform an expansion around this constant value. 
Keeping the first non-trivial term, 2(r, a) can be approximated by the following expression: 

Zir, a)^Z,- tl * (2.20) 
while for the dilaton $ we have: 

p3 ^3/2 

^-mr,a) ^ ^ ^* (2.21) 

8 [Vr^Z. + a^j' 
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Notice that the expression for Z in fl2.20p gives the exact result (12.161) for cr = 0. The above 
analysis suggests that in the UV region the combination r^Z^ + plays a significant role. 
Having this in mind we define a new set of variables, u and d, as follows: 



u 



a/ r^Z* + cr^ & tan a 



a 



'Z^r 



with 



TT 

< a < - . 



f2.22) 



The function Z as well as the dilaton generally depend on both coordinates, u and d, but 
in the UV limit of large u the a dependence disappears. Actually, their expressions when 
M — >• cxD are: 



z ^ z, & 



-2* 



2v^ \ u) 



(2.23) 



Using these values in the metric ansatz we end up with the following expression: 



Ubjjy 



2^2 /M\3/2 

^ \r) 



c?xi_2 + Z^R^dQl 



+ 



2V2Z. 



1/4 



-) 



3/2 



1/4 



u 



1/2 



(id^ + cos^ a d^ll + sin^ d (dip + cos ^ c/0)" 



(2.24) 



where dQ"^ 



dO"^ + sm^ 9 dcj? is the line element of the (^, 0) two-sphere. In order to interpret 
the meaning of the results just found, let us recall that, given a background of type IIA 
theory such as the one we are considering, one can generate a solution of eleven- dimensional 
supergravity by uplifting the metric by means of the standard formula: 



ds\^ 



'dsl, + ea* 



(dz) 



(2.25) 



where z is the eleventh M-theory coordinate. We shall apply (I2.25P to the ten-dimensional 
UV metric and dilaton written in eqs. (12.241) and (12.231) . After changing the radial variable 
M by a new coordinate y, defined as: 

2R 



y 



z^ 



u 



(2.26) 



the resulting eleven-dimensional UV metric becomes: 

,2 



ds^ 



11 



i?2 



"2 3 + dnl 



dy-^ 2 

y 



+ 



i?2 



dce'^ + cos^ d cifij + sin^ d (dip + cos 9 dcj))' 



(2.27) 



where dx\j^ = dx\2 + dz"^. From (12.271) we conclude that the uplifted metric is of the form 
AdS-j X S^, with the AdS^ having two of its directions compactified in a two-sphere and with 
the S"^ being fibered over this S"^, Notice also that the radius of the AdSi is just i?, whereas 
the has radius R/2. These results are, of course, consistent with the origin of the solution 
as coming from M5 wrapped on a two-cycle. 
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3 Addition of flavor branes 



In this section we will start exploring the possibility of finding the dual of the A/" = 4 3d 
gauge theory with matter hypermultiplets in the fundamental representation of the gauge 
group. We will achieve this by adding fiavor branes to the setup of section [2l These fiavor 
branes will be D4-branes extended along the three Minkowski directions and as 

well as the ip and a directions of the Calabi-Yau. At the same time they will be located 
at particular fixed values of the S"^ sphere and of the transverse M'^, as represented in the 
following array: 
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N2 




D4 (color) 
























Nf D4 (fiavor) 























It is worth pointing out that, in the above setup, the fiavor branes wrap a non-compact 
direction of the internal Calabi-Yau. This is, actually, a standard requirement which one 
should demand of these kinds of setup in order to convert, in the appropriate decoupling 
limit, the gauge symmetry living on the worldvolume of the fiavor brane into a global (fiavor) 
symmetry. 

To describe more precisely the embedding of the fiavor brane let us choose the following 
system of worldvolume coordinates: 

r = (x°,x\a;2,a,^) . (3.1) 

Then, the embedding is determined by the condition that all other remaining ten-dimensional 
coordinates are constant. By analyzing the kappa symmetry of the worldvolume action of 
the fiavor brane we will explicitly check in appendix [B] that the embedding just described 
preserves the same supersymmetry as the unfiavored setup. Moreover, the position in r 
of the fiavor brane, rg, represents the distance between the two sets of branes and has a 
well-defined meaning in the gauge theory dual. Indeed, tq is related to the mass uiq of the 
matter hypermultiplet by means of the following relation: 

(3.2) 



2'Ka' 



In the so-called quenched approximation the effects of the quark loops in the field theory 
observables are neglected. This approximation is well justified when Nf <^ A^^^ and cor- 
responds, on the gravity side, to treating the fiavor branes as probes and neglecting their 
infiuence on the metric. Later in this paper we will make use of this approximation to study 
several aspects of the gauge theory, such as the meson spectrum and the quark-antiquark 
potentials. However, in the remainder of this section we will analyze, in a certain approxi- 
mation, how the backreaction of the fiavor branes modifies the solution described in section 

El 
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3.1 Including the backreaction 



Let us study the backreaction of the flavor branes on the background in the case in which 
the number of flavors Nf is large and of the same order as the number of colors Nc- From the 
field theory point of view this limit was considered a long time ago by Veneziano [6]. Here 
we will follow the approach pioneered in ref. [Tl], which is based on the observation that, 
when Nf oo, one can homogeneously distribute the Nf flavor branes in their transverse 
directions (for a clear discussion on the validity of the DBI+WZ action for a large number 
of smeared branes, see section 7 of [21]). Notice that, when the branes are embedded 
as explained around (13. ip . they preserve the same supersymmetries independent of their 
position in the transverse space. Actually, we will consider a distribution of branes with 
a fixed value rg of the r coordinate and smeared along the angular coordinates {9, 0) and 
{6, (f). In order to figure out how this smearing is implemented, let us recall that the action 
for a stack of Nf D4-branes is given by the sum of the DBI and WZ terms: 



S 



flavor 



[ c?'ee-* V-detGs + T,J2 I ^' 

J Mr, J Mr, 



(3.3) 



with Gs being the induced metric on the worldvolume and 6*5 the pullback of the RR four- 
form potential to A^s. The smearing procedure amounts to promoting the infinite sum 
appearing in the action to a ten-dimensional integral. For the WZ part we have: 



E 

Nf 



Mr, 



QAC. 



(3.4) 



Mio 



where 17 is a five-form proportional to the volume form of the transverse space, namely: 



N. 



IGtt 



— 6{r - tq) dr a UJ2 a U2 



(3.5) 



with UJ2 being given by (12. 5p and ijJ2 = sin 6d6 A d(j). The smearing form Q is normalized as: 



Q = Nf. 



(3.6) 



Similarly, to obtain the smeared version of the DBI part of the flavor brane action we perform 
the substitution: 



J2 [ rf^e e"* \/- det 65 ^ / d^°x e"* V- det G | | , (3.7) 



where is the modulus of fl: 



\ 



5! 



s n 



k=l 



(3.8) 
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Therefore, the smeared DBI+WZ action of the flavor branes is: 



S 



flavor 



-T4 f d^^x e"* V- det G \n \ + f ^] A C5 . 



(3.9) 



By inspecting the WZ term in the action (13. 9p one readily concludes that the flavor brane 
acts as a source for the RR six-form Fq = dC^ which, in turn, induces a violation of the 
Bianchi identity of F4 = *Fq. Actually, one can prove that this modified Bianchi identity 
becomes: 



dF^ 



2Nr 8 



5{r — tq) dr /\ UJ2 /\ UJ2 ■ 



(3.10) 



where, in the last step, we have used that 2k^q = {2T[y gl{a'Y and T4 = 1/ {2TiYgs{a'Y^ 
and, thus, 2k\q = n"^ R^/Nc- 

Let us now formulate a new ansatz for this backreacted flavored setup. First of all, we will 
adopt the same ansatz (12. ip for the ten-dimensional metric. However, as is clear from (I3.10p . 
we should change the ansatz for F4 in order to reproduce the modified Bianchi identity. 
Actually, the natural ansatz both satisfying (13.101) and generalizing (12.60 is: 

~ ~ ( Nf \ ^ 

F4 = -{gda + g'dr) A 002 A {d'^ + cose d(f))+ lg + ~ ^^^^2, (3.11) 

where is the Heaviside step function. Proceeding as in the unflavored case and substituting 
the new ansatz (13.111) for F4 into the equations for the supersymmetric variations of dilatino 
and gravitino we have: 



9 + 



2N^ 8 



e(r - tq) 



- R^r'^Z' 



-4$ 



a = R^ZZ 



(3.12) 
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-o 



R-^Z-^e-^" 



9 + 



Nf R' 



2Nr 8 



e(r - rq) 



The set of projections on the Killing spinors needed to arrive at (I3.12p is the same as in the 
unflavored case. Thus, the flavored solutions preserve also the same eight supersymmetries 
as the unflavored ones. Moreover, as in (12. 8p . the last equation in the system (13.120 is not 
independent of the others. One can also verify (see appendix |A]) that the equations of motion 
of F4 and the metric (including the contribution of the DBI action to Einstein's equations) 
are satisfied if the system (I3.12p holds. The analogue of the PDE (12. 9p for the flavored case 
is: 

'2 ' rZ" + 2Z') + ^^-^^^ 



r'Z [Z - oZ 



r a { r Z 
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As in the unflavored case, we can integrate the function Z for cr = 0. Indeed, it follows from 
f l3.12p that g is independent of r when the variable a vanishes. If we call (70 this constant 
value of g then, for r > tq, we have: 



Qo H 



1 



and after using for go the same value as in the unflavored case, we get: 



Z'{r,0) 



R 



2N, 



r > 



(3.14) 



(3.15) 



When r < tq, by simply putting Nf = on the right-hand side of fl3.15p . we recover the 
unflavored result. Integrating fl3.15p and imposing continuity for the solution along r = tq 
we have: 



Z(r,0) = 



R Nf R 
— Q[r — rn) 



2 Nr. 



e(r-rQ) 



(3.16) 



where the constant Z^, is the same as in (12.160 . We have not been able to obtain analytically 




Figure 1: Solutions for Z for different values of X — as (X function of r and a. These 
curves correspond to k = Z^ = 1 and tq = 5. 



the solution of the flavored system fl3.12p for arbitrary values of r and cr. Instead, one 
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can integrate numerically this system of equations. Since the system (13.121) reduces to the 
unflavored one written in (12.81) . we can assume that the solution of (I3.12p reduces to the one 
written in (I2.12I) - (I2.14I) for r < rg, while at r = the functions 2, g and $ are continuous 
and Z' has a discontinuity independent of a and given by: 

Z'ivQ + e, a) - Z'ivQ " ^) = ^ ■ (3-17) 

Using the continuity and change in derivative of Z at r = rg and 13.161 as the boundary 
conditions, we are able to numerically integrate the equation of motion for Z up to a finite 
value of r after which the solution becomes highly oscillatory and unstable. The result of the 
numerical integration for Z{r, a) is shown in figure [3?T] for several values of N^. and Nf. It is 
quite evident from figure 13.11 that Z has a wedge shape at the position r = rg of the flavor 
branes. This means that there will be a curvature singularity at this point, which is actually 
needed to match a similar term in the Einstein equations coming from the energy-momentum 
tensor of the flavor brane sources. As in the unflavored case, the geometry is also singular 
in the IR due to the enhangon phenomenon. Outside these regions the metric is regular 
and one can verify that the components of the Ricci tensor depend on N^. and Nf through 
the combination Nc — ^f, which scales as N^. in the Veneziano limit N^Nf — > oo with 
Nc/Nf fixed. Thus, outside the location of the source, the region in which the supergravity 
approximation is valid is not modified by the backreaction of the flavor branes, as happened 
in the backgrounds of refs. [TT], [13 [22] ■ 



4 Connection with gauge theory 

The Af = 4 three-dimensional gauge theories can be regarded as the reduction to three 
dimensions of the M = 2 gauge theories in four dimensions [27]. The field content of the 
different 3d supermultiplets can be easily obtained by dimensional reduction of the one 
corresponding to A/" = 2, ci = 4. Let us consider, for example, the vector multiplet which, 
in four dimensions is composed by a vector field, a Dirac spinor and a complex scalar. By 
reducing to three dimensions one gets one vector, two Dirac fermions and three scalars. 
Similarly, one can verify that the N' = 2,d = 3 hypermultiplet contains two Dirac fermions 
and four complex scalars. Given this field content it is straightforward to find the one 
loop running coupling constant gvuilj) of the theory (see, for example, the appendix C of 
[lOj). Indeed, let us consider a SU{Nc) gauge theory with Nf matter hypermultiplets in the 
fundamental representation. Then, one can show that: 



1 



9ym 9ym 



9YMNcf, Nf 



(4.1) 



AiTfi V 2Nc 

where fi is the energy scale. The one-loop result displayed in eq. (14. ip is, actually, exact 
in perturbation theory. Notice that (14. ip shows that the Af = 4 theory has negative beta 
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function when Nf < 2Nc, while for Nf > 2Nc the beta function changes its sign and becomes 
positive. In the borderhne case Nf = 2Nc the one-loop beta function vanishes and the 
coupling does not run anymore in perturbation theory. In the next subsection we will show 
that our gravity solutions match perfectly the behavior f l4.ip . both in the unflavored and 
backreacted flavored cases. 

Besides the perturbative result just reviewed, the A/" = 4 3d theories have a very rich 
non-perturbative structure. Indeed, the Coulomb branch of vacua of these theories is a hy- 
perkahler manifold, which is isomorphic to the moduli space of three-dimensional monopole 
solutions of a different gauge theory [27ll2H]- Moreover, the Higgs branch is also a hyperkahler 
manifold. Furthermore, these theories display the phenomenon of mirror symmetry [29| [30], 
which is a duality between two different AT = 4 3d gauge theories which exchanges the Higgs 
and Coulomb branches, as well as the Fayet-Iliopoulos and mass terms. When these A/" = 4 
theories are realized in the type IIB string theory as in [28], the mirror symmetry is just a 
manifestation of the underlying S-duality of the type IIB string theory (see ref. [31] for a 
review) . 

4.1 Probe calculation 

In order to extract information about the gauge theory living on the D4-branes, we will 
study, following ref. [10], the dynamics of a color D4-brane probe wrapping an S"^ and 
moving under the influence of the metric and RR form of the background. The action of 
such a probe will be: 

S = -Ti j rf^^e"* y^-det(G'5 + 2na' F) + Ta j (c^ + 2^(73 A , (4.2) 

where (a = 0, ■ ■ ■ , 4) is the set of worldvolume coordinates along which the color D4-brane 
is extended, F is the field strength for the worldvolume gauge field and the hat over G5, 
C5 and C3 denotes the puUback over the worldvolume of the D4-brane. In particular, let us 
evaluate the action (14.21) in the case in which the ten-dimensional metric is of the form fl2.ip . 
Let us choose the worldvolume coordinates as = (x°, x\ x^, 6^, 0) and let us assume that 
the brane in embedded in such a way that the remaining ten-dimensional coordinates are 
constant and the worldvolume gauge field F is zero. Then, the determinant of the induced 
metric for such a configuration takes the form: 



e-^^-deiG, = R'Ze'" sin0~ ^1 + cr^ ■ (4-3) 

For the configuration we are considering, the only non-vanishing contribution to the WZ 
part in (14.21) is the term containing C5. The expression of this five-form potential for our 
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solutions has been evaluated in appendix |X] (eq. (1A.9P ). The corresponding puUback to the 
D4-brane worldvolume is: 



Cr> = R^Z e^* sin 9 dx° A dx^ A dx^ A dO A d4) . (4.4) 

Let us now substitute (14.31) and (14. 4p into (14.21) . The result is just minus the static potential 
between the stack of A^^^ color branes and the probe, namely: 



^ , = -T4 / d^xdMd) R'^Ze^'^ sine \ 1 + a 



cot^ I 



pot — —^4 / u. j^uuuyj iL ^ c oiii u |^ y ± -r u ^2 ^2 ^4<i> 



(4.5) 



Notice that the right-hand side of (14.51) only vanishes when cr = 0, which should be inter- 
preted as the point of the Calabi-Yau in which one can place a color brane without breaking 
super symmetry. Actually, one can check this statement directly by studying the implemen- 
tation of kappa symmetry for the different D4-brane embeddings (see subsection 14.21 and 
appendix [B| . Let us thus assume that our probe brane is located at a = and that we 
switch on a worldvolume gauge field F^iy whose only non-vanishing components are those 
along the Minkowski directions x^. We will expand the DBI lagrangian density (integrated 
over the angular directions 9 and 0) up to quadratic order in the gauge field F^^. The result 
of this expansion can be parameterized as: 

[ ded4>CDBi = --^i— tr[F^,Fn , (4.6) 

J quadratic ZgyMyl'^) 

where gvui^J') is, by definition, the Yang-Mills coupling at the renormalization scale /i. 
Actually, by performing explicitly the calculation, one gets: 

^ Z(r,(T = 0). (4.7) 



Clearly, the variable r should be related to the energy scale n of the field theory. The natural 
radius-energy relation is given by: 

r = 27ra'/i, (4.8) 

which we will assume to hold in the following. Let us now suppose that r > tq or, equiv- 
alently, that /i > mg, where ttiq has been defined in (13.21) . After substituting in (14. 7p the 
value of Z{r,a = 0) as given by (13.160 for the general flavored solution we have: 

f. ^ (4.9) 



Moreover, by defining the bare UV YM coupling as: 



1 Z, i?2 N 



9ym 2 71 gs{a'y/^ SirrriQ 



(4.10) 
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we get that eq. (14. 9 p matches perfectly the field theory expression (14 .11) . Moreover, by 
allowing the brane probe to move in the transverse fiat space, and by looking at the action 
of the transverse scalar fields obtained by expanding the DBl+WZ action, one can obtain the 
metric of the moduli space in the Coulomb branch. After dualizing the worldvolume gauge 
field one can check that the moduli space is indeed hyperkahler and, actually, its metric has 
the form of the Taub-NUT space [10]. 

4.2 Higgs branch 

The Higgs branch of the A/" = 4 gauge theory is a phase in which the quark hypermultiplets 
acquire a non- vanishing expectation value. On the field theory side one can study the theory 
in the Higgs branch by turning on an extra Fayet-llioupoulos coupling in the lagrangian. 
The realization of this mechanism in a brane setup is well-known [28] (see [31] for a review). 
Indeed, as argued in ref. [28], one should reconnect the color and fiavor branes in a super- 
symmetric way. In our holographic setup we should look for D4-brane embeddings which 
are compatible with all the supersymmetries of the gravity solution and such that they can 
be interpreted as representing a recombination of color and fiavor branes (see [22l [32| l33] for 
a similar analysis in other brane setups). Recall that both types of D4-branes are extended 
along different directions of the Calabi-Yau cone. Indeed, the color branes are extended 
along (^,0) at o" = while the flavor branes extend along {(T,tp) at fixed angles {0,(f)). In 
order to find a configuration interpolating between these two situations it is natural to use 
the same system of worldvolume coordinates as in (13.11) and look for an embedding such that 
6 and are no longer constant. Instead, they will depend on the other coordinates of the 
CY2, namely: 



In order to determine the D4-brane embeddings of the form (14. lip which preserve the su- 
persymmetries of the background, one has to study the kappa symmetry of the brane probe. 
This analysis is performed in detail in appendix [Bl The final result found in this appendix 
can be nicely recast in terms of the following two complex coordinates of the CY2: 



It turns out that any holomorphic embedding of the type Ci = f{.C2), with / arbitrary, solves 
the kappa symmetry equations and, thus, preserves the supersymmetry of the background. In 
order to make contact with the field theory analysis, it is rather natural to restrict ourselves 
to embeddings characterized by a polynomial equation of the type: 

Cr = C , (4.13) 

with C being a (complex) constant and the exponents pi and p2 are constant integers. Notice 
that, in the Coulomb branch, the embedding of the color branes correspond to eq. (14.130 



= 0(f^,^) • 



(4.11) 





(4.12) 
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with pi = 0, p2 = 1 (and C = 0), while that of the flavor brane is given by pi = 1, p2 = 0. It 
is thus natural to think that the Higgs branch embedding we are seeking is the one obtained 
by taking pi = 1, p2 = I in (14.131) . Notice that, in this case, the functions ()4.11l) are just 
given by: 



2 ^* 



sin^y = ^, = V + Vo, (4.14) 

where and ipQ are constants (related to C in (14.131) ). Notice that a^, is just the minimal 
value of the coordinate a (which occurs for 6 = n) which should correspond, in the field 
theory side, to the Higgs VEV. 



5 Meson spectrum 

In this section we analyze the mass spectrum of the mesonic excitations for the backgrounds 
of sections [2] and [3l In order to address this question we will add a flavor D4-brane probe to 
these backgrounds and we will study the normalizable fluctuations of its worldvolume fields. 
This analysis will lead us to a problem of the Sturm-Liouville type with a discrete set of 
eigenfunctions and eigenvalues. 

Let us consider a flavor D4-brane in a geometry of the form (12. ip . We will choose the 
system (13. ip of worldvolume coordinates and we will consider an embedding in which the 
radial coordinate r and the angular coordinates 6, (p, 9 and are constants. Let us denote by 
Tq the constant value of r for this configuration (r^ is related to the mass of the external 
quarks as in ( 13. 2p . i.e. rriq = rq/2iTa'). The induced metric for such a configuration is: 

= e'*(''"") rfx?2 + -T7 r da^ + a^d^? . (5.1) 

Z{rg,a) L J 

Let us now perturb this static configuration by deforming it in the radial direction as: 

r = Tg + f{x'^,a,ilj) , (5.2) 

where the fluctuation f is small. In all the calculations of this section we will take Vg such 
that the r = surface does not enter the enhangon region and, thus, a = is the minimal 
value of the coordinate a in the worldvolume. This can always be achieved by taking Vg 
sufficiently large. 

Of course, the perturbation (15.21) is not the most general one. However, one can check 
that, at quadratic order, f does not mix with other fluctuations and, therefore, it can be 
studied separately. By expanding the corresponding DBI+WZ action one can verify that. 



^One should take into account that, due to the different signs in the two projections in (jA.4p . the WZ 
coupling of the RR potential C5 for the flavor brane must be opposite to the one appearing in 
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up to quadratic terms, one gets that the lagrangian density for f is: 



4 



{d^,rf + a{d^rf + -{d^rf , (5.3) 



e4*Z ■ ■ -a 



where $ and ^ should be understood as functions of a at r = r^. The equation derived from 
the lagrangian (15.3P is given by: 

9^ H.r] + ^l^f- + -^l^ = • (5-4) 

e is cr ^ 

To find the solutions of this equation, let us separate variables as: 

f = e^'^^ e^''^ , (5.5) 

where I is an integer (which can be taken to be non-negative without loss of generality) and 
is a momentum along the Minkowski directions x^. Let us also define as = — fc^, 
where the square is computed with the fiat Minkowski metric in 2 + 1 dimensions. Plugging 
(15. 5p into (15. 4p . we arrive at the following equation for xi'^)'- 

[ad^x] + \^ - -jx = , (5.6) 

Interestingly, by means of a suitable change of variables, the fluctuation equation (15.61) can 
be written as a Schrodinger equation. Indeed, let us define the variable y as follows: 

= a . (5.7) 

Notice that y G (— oo, +oo). In terms of ?/, the equation ( 15. 6p of the fluctuations can be 
written as: 

0-ni/)x = O, (5.8) 
where the potential Viy) is given by: 

V{y) = ^"-M'^. (5.9) 

From the reformulation of eqs. (15.81) and (15.91) of the fluctuation equation, one can easily 
obtain the asymptotic value of x when y ±oo. Indeed, from the behavior of $(rg, a) 
and Z{rg, a) when o" — > 0, oo one easily gets that V P when y —>■ ±oo. It follows that, 
in these asymptotic regions, the two independent solutions of (15.81) when / ^ are just 
X ~ e^'^ = 0"^' (when / = these solutions behave as x = constant, log cr when a ^ 0, oo). 
The normalizable solutions, which can be identified with mesonic excitations, are those that 
are regular as a = 0, oo. These solutions only exist for some discrete set of values of the mass 
M, which can be determined numerically by means of the shooting technique (see below). 
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The previous analysis only applies to the particular fluctuation f in the radial direction 
of the transverse M^. However, one can check that the same equation (15.61) describes the 
fluctuations of the other two coordinates 9 and (j). Furthermore, in order to study the 
fluctuations of the coordinates 9 and that determine the position of the flavor brane in the 
CI2, let us write: 

^ = ^0 + ^(a;^ ^, a) , = 00 + , (5.10) 

where ^0 and 0o are the constant unperturbed values of the angles 9 and (p. By plugging (15.101) 
into the DBI+WZ action and expanding the result up to quadratic order in the fluctuations 
one gets a lagrangian in which 9 and are coupled. It turns out that, as happened in refs. 
p!5] and [22], they can be easily decoupled by defining new fields x+{^) and X-{^) as follows: 

^ = i (x+(a) + x-{<y)) sin 00 e^'^ sin(/^/.) , 

- e^'" cos(Z^) , (5.11) 

where / is a non-negative integer. One can verify that, indeed, and satisfy two 

different decoupled equations. Actually, it is possible to perform a further redefinition to 
the functions x± such a way that the new functions satisfy the scalar fluctuation equation 
(15. 6p . Similarly, one can show, as in ref. [15], that the fluctuations of the worldvolume gauge 
field also satisfy (15. 6p . These facts imply that, up to finite shifts in the quantum numbers, all 
fluctuation equations lead to the same set of eigenvalues, as expected from supersymmetry. 
Therefore, one can concentrate on studying the basic equation (15. 6p . This is what we will 
do in the next two subsections. 



5.1 Quenched mesons in the Coulomb branch 

Let us specialize the previous analysis to the case in which the background geometry is the 
unflavored one of section [2] and we add a flavor D4-brane probe. In this case the function 
Z and the dilaton $ are given by eqs. (12.131) and (I2.14p respectively. The numerical results 
obtained by means of the shooting technique are shown in figure |2l In general, for given 
and /, one obtains a tower of discrete normalizable states labelled by a principal quantum 
number n. It is interesting to notice that, for fixed n and Z, the masses grow with as 

One can get a rather accurate estimate of the mass spectrum by applying the semiclassical 
WKB method to the Schrodinger problem (15. 8p (see refs. [3l]-[37]). By applying this method 
to our particular case, we get the following mass formula: 

= -f-T2 + + 3/) , (5.12) 
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Figure 2: Mass eigenvalues from the improved WKB formula (eq. 15.16)) (dots) and from the 
numerical calculation (solid line). The numerical calculation has been performed with the 
unflavored background with k = = 1. 



where n is a non-negative integer > when / ^ and > 1 for / = 0). In eq. (15.121) 
C{fq) is the following integral: 



+00 



da 



^g4*{r„<x) z{rq, a) 



(5.13) 



Notice that C(^(?) determines the mass gap of the mesonic spectrum. One can get an approx- 
imate expression for C(^g) by using inside the integral on the right-hand side of (15.131) our 
asymptotic expressions (I2.20p and (I2.2ip . By doing so one arrives at the following analytic 
expression for C('^(j): 




i?3 1 



+ 00 



da 



(^2 + Z^rD'i 



7r3 i?3 
16 



r(!) 



(5.14) 
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It is interesting to stress that, within this approximation, ({rg) is independent of 2^*. Ac- 
tually, by using the value of R written in (12.21) . as well as the relation between and the 
quark mass rrig (see (13. 2p ). one can write the following WKB formula for the meson masses: 



r(!)l' 



WKB 



TcgsNcVa' 



n+1 



n 



31 



(5.15) 



This formula reproduces rather well the numerical values of the masses for / = 0, 1 and, 
actually, captures accurately the meson mass gap and its dependence on the quark mass rriq. 
However, the degeneracies observed in the numerical results when one varies the quantum 
numbers n and I are not reproduced by fl5.15p . Nevertheless we found that the numerical 
results are fully recovered if (15.151) is changed to: 



r( 



1 4 



TrgsNc 



'a' 



n + l] {n + 2l + l 



(5.16) 



The comparison between the predictions of (15.161) and the numerical values of the masses is 
shown in figure [2l It follows from the inspection of this figure that (I5.16P provides a very 
good fit of the masses found by the shooting technique. 



5.2 Unquenched mesons in the Coulomb branch 

We can perform the analysis of the meson spectrum by looking at the fluctuations of a flavor 
brane probe in the background given by the backreacted flavor branes. This is equivalent 
to studying the fluctuations of a single quenched flavor in the presence of Nf unquenched 
flavors. In the large Nf limit this should be an accurate approximation. What is interesting 
to note here is that because the functions Z and $ are themselves continuous in r and are 
being taken at the value r = r^, the eigenvalues for a flavor probe at this position will be 
exactly the same as that for a probe without the backreacted flavor. 

Thus, by inspection we see from the form of the equations of motion for the probe brane 
fluctuations that the spectrum will be identical in the flavored and unflavored case if the 
probe brane is placed at the same position as the backreacting branes. The difference between 
the quenched and unquenched case will only be felt if we introduce a flnite distance between 
the probe and backreacting branes, corresponding to Higgsing a U{1) of the U{Nf) flavor 
group. 

The case where the probe and backreacting branes are separated is considerably more 
complicated, computationally, because in order to calculate the meson spectrum it is vital 
to have a numerically stable geometry. The calculation of the function Z from the partial 
differential equation is inherently an unstable calculation and therefore the calculation of 
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Figure 3: Variation of the meson spectrum calculated on a probe brane at Vg = 5 with 
backreacting flavor branes at varying positions tq. The four lines on each plot correspond 
to Nf/Nc = 1 (flattest line), Nf/N^ = 2, Nf/N^ = 5 and Nj/N^ = 10 (least flat line). It is 
clear to see that each of these return to the quenched case in the case where the backreacting 
branes lie close to the probe brane. All calculations have been done for k, = = 1. 



the spectrum is difficult. We have however managed to calculate the spectrum in a narrow 
range of parameter space where the probe and backreacting branes are not very far apart. 

In figure [3] we study the effects of placing the backreacting branes at a variable position, 
given by rg while always keeping the position of the probe brane at = 5. On each graph 
the lines correspond to Nf/N^ = 10, Nj/N^ = 5,Nf/Nc = 2 and Nf /N^ = 1. Clearly the 
Nf/Nc = 1 line is the flattest in each case and the Nj/Nc = 10 line has the most variation, 
as expected as the geometry is most altered with a higher ratio of flavor branes to color 
branes. We have checked that the curves in figure [3] can be fitted to an expression of the 
type = a + Nfb [y- — l) , where a and b are coefficients that are independent of Nf. 

5.3 Mesons in the Higgs branch 

We will now consider the fluctuations around a non-trivial embedding of the type studied in 
subsection 14.21 (eq. ( 14.13^ ). Recall that we argued in this subsection that these embeddings 
correspond to configuration in which both types of D4-branes (color and flavor) are recom- 
bined, realizing the Higgs branch in our brane setup. More concretely, we will concentrate 
on studying the embedding (14.13^ for pi = p2 = 1, which can be rewritten as in (I4.14p . We 
will continue to use fl3.ll) as our system of worldvolume coordinates but now 6 and will 
no longer be constant but given by the functions displayed in (14.141) . For simplicity we will 
restrict ourselves to studying the fluctuations of the radial coordinate r around a fixed value 
Tq (see eq. (15. 2p ). After separating variables as in (15. 5p the fluctuations equation can be 
written in the Schrodinger form (15.81) . if we introduce a new variable y as: 

= 0" — cr* , —oo < y < +oo , (5-17) 
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where a^, is the minimal value of the coordinate a for the embedding (see eq. fl4.14p ). Recall 
that cr* parameterizes the Higgs VEV. The Schrodinger potential V is now given by: 



V = - - M'' 
4 



e^y o2 Z 



(5.18) 



Note here that in the limit / 21 and — we recover the Coulomb branch potential. 
However, one can check that this limit is not continuous and the last term in (15.181) means 
that the Coulomb spectrum is never recovered taking this limit in the Higgs branch. 

In order to have an idea of the mass spectrum associated to the potential (I5.18p . let us 
apply the WKB method. The mass gap is now determined by an integral similar to (I5.13p . 
which is now given by: 




When cr=K is not very small the second term inside the square root in (15.190 dominates the 
integral. Actually, by applying the same approximations that led to (I5.14p . one arrives at 
the following approximate expression for ({rq,a^:): 



C(rg,a,) ^ -kR^Z, . (5.20) 

Notice that C in (]5.20p is independent of and a* and only depends on the asymptotic value 
of Z^,. Actually, after improving the WKB formula as in (I5.16p . we arrive at the approximate 
formula for the mass levels in the Higgs branch: 

^maas ^ {n + l){n + l + l) , (5.21) 

which reproduces rather well the numerical values. 



6 Wilson loops and the qq potential 

In order to investigate the unflavored background more, we now turn to another direction 
of a more phenomenological nature, the Wilson loop operator. On the gravity side, the 
Wilson loop expectation value is calculated by minimizing the Nambu-Goto action for a 
fundamental string stretching into the dual supergravity background, whose endpoints are 
constrained to lie on the two sides of the Wilson loop. Below, we first briefly review the 
procedure for calculating Wilson loops in supergravity. We will discover that the qualitative 
behavior of this observable critically depends on whether the constant k, appearing in the 
function (12.120 . is smaller or larger than 1/16. We will discuss these two cases separately in 
two different subsections. 
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6.1 General formalism 



As stated above, the calculation of a Wilson loop in the gravity approach amounts to extrem- 
izing the Nambu-Goto action for a string propagating in the dual geometry whose endpoints 
trace the loop. To describe such configuration in our setup let us choose the time t and a 
Minkowski coordinate x as worldvolume coordinates of the string and let us consider the 
trajectory with: 

r = r{x) , a = , rest = constant . (6-1) 
Then, the induced metric on the string worldvolume can be easily found from (12.11) . namely: 

e^-f M) [ _ ^ ^ g-4<i.M) ^/2) ^^2 ] ^ 2) 

with $(r, 0) being given by (I2.17p . which we can rewrite in terms of the enhangon radius 
defined in (EH]) as: 



^2 / ^2 



2rrp + 16 /tr? 



The Nambu-Goto action for the configuration (16. ip is: 

S = — [ dtdx Ve^*(^'°) + r'2, (6.4) 
27r J 

where we have taken a' = 1 and the prime denotes a derivative with respect to x. Since 
the action does not explicitly depend on x, the system has a first integral tq, which can 
be identified with the turning point of the solution. Solving the corresponding first-order 
equation for x in terms of r we find that the linear separation of the quark and antiquark is: 

roo g2<J>{ro,0) 

L = 2 dr . (6.5) 



ro e 



(r,0) ^g44.{r,0) _ g4$(ro,0) 



The energy of the configuration can be obtained from the action (16.40 . Subtracting the 
self-energy contribution, we obtain 

I r<x> g2*(ro,0) ]^ roo 

E = - dr , - - / rfr , (6.6) 

TT^ ^g4<I.(r,0) _ g4#(ro,0) yr ^^^^^ 

where rmin is the minimum value of r allowed by the geometry. In the specific cases, we are 
supposed to solve for the auxiliary parameter tq in terms of the separation length L. Since 
this cannot be done explicitly in practice, one regards (16. 5p as a parametric equation for L in 
terms of the integration constant rg. Combining it with (16. 6p for E, one can then determine 
the behavior of the potential energy of the configuration in terms of the quark-antiquark 
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Figure 4: Effective tension of the string as a function of r. On the left is the case with 
/t < 1/16 and on the right with n > 1/16. 



separation. Using fl6.5l) and fl6.6l) we arrive at the following relation which will prove useful 
in the following: 

^ = ■ (6-7) 

The minimum value of r allowed by the geometry depends on k. Indeed, as argued at the 
end of subsection 12.11 the numerator in the right-hand side of (16.31) has real roots only when 
K < 1/16. In this case r^m = th, where th has been written in (12.191) . Conversely, when 
K > 1/16 the minimum value of r is = ^e- 

Following the analysis in [21] , we can determine the qualitative behavior of the system sim- 
ply by looking on the asymptotic expansion of the dilaton in (16.31) . Utilizing their formalism 
we write the lagrangian as: 

C = T{r)^/1 + Y{r) r'\ (6.8) 

where T(r) is the effective tension of the brane depending on the dilaton. In our background 
this is given by: 

T(r) = ^ = e-(-°) . (6.9) 

When K < 1/16 both the tension of the string and Y behave in the IR with a power law, 
T = yT72 ~ (r — rminY^'^- Following the reasoning of ref. [21], due to the value of the 
exponent of Y, the string will not reach an infinite length in the x direction when it explores 
the far IR of the dual theory. As we will see in the next subsection this is verified both by 
our numerical and analytical results. 

When K > 1/16 the string feels an effective IR wall at the point r = r^, where the tension 
is minimized, so T'{r^:) = and T"{r^) > (see figure H]). The actual value of can be 
obtained analytically by finding the roots of a cubic polynomial. These roots depend on 
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Figure 5: k < 1/16: Separation length between a quark and a antiquark as a function of tq 
and energy as a function of the separation length. The dashed (solid) curve is the case for 
K = 1/16 {k = 1/18). 

K and Te and one can check that r^, 2re when k is large. Since the background at the 
point r = is smooth and the functions of the metric do not vanish, they admit a Taylor 
expansion as: 

T = n + T2ir-ny + ..., T2>0, 

Y = Y, + Yi{r + (6.10) 

Following again the reasoning in [21], when the tip of the string comes close to the IR wall 
the length will diverge. Again this is something we will verify in the next subsection. 

6.2 Case with < 1/16 

For large values of tq the potential is, as usual, Coulombic while their exists a maximal 
separation k < 1/16. At the values of tq where the separation is maximal, the energy is also 
maximal and always positive (see figure [5]) . This means that there is a screening behavior 
because as the potential turns positive, a configuration of two separate strings is energetically 
favored and corresponds to a vanishing force between the charges. As ro Tmin = i^h we 
find by studying the behavior of the integrals in eqs. (16. 5p and (16.61) in this case that both 
length and energy approach zero as: 

L^V2R f^V'\/^ In 
\ rn J \ r-H 

and 

E ^ In 

27r 



To - Th 



(6.11) 



(6.12) 
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Figure 6: k > 1/16: Separation length between a quark and an antiquark as a function of 
tq and energy as a function of the separation length. The left (right) curves correspond to 
the case in which k is below (above) the critical value and the central curves are the case in 
which K, has exactly the critical value, Hcr — 63.8. 



6.3 Case with k > 1/16 



In this case we follow the same steps as before but now we should have in mind that 
rmin = As usual, the dependence of the potential energy for small separations L of 
the quark-antiquark (corresponding to large tq) is Coulombic. In order to have an approx- 
imate expression for the separation length and the energy we have to move to the region 
of the parameter space where k is large. In this case the effective IR wall is located at 
r = r^: ^ 2re and the potential gives a linear confining behavior with: 



1/2 



V3 + 4k 
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6r^ 



ro - 2re 



and 
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4r^ 
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3 + 4k 
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6r. 



K 




(6.13) 



(6.14) 



ro - 2re_ 

For intermediate values of the separation length, the behavior of the system depends crucially 
on the value of n. There is a critical value Kcr — 63.8, such that for k, > the behavior of 
the length and energy curves resemble the Van der Waals isotherms for a statistical system, 
with ro, L and E corresponding to volume, pressure and Gibbs potential respectively (see, for 
instance, [3H]). In the following we will describe the situation referring to the corresponding 
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plots in figure [6l In the region below the critical point, k < Kct, the energy is a single- valued 
function of the length, connecting a Coulombic with a linear confining phase. For values 
above the critical point, k, > Hcr, the Coulombic phase for the energy at small distances is 
followed by a situation where the energy becomes a triple- valued function of the length, with 
the physical state being the one of lowest energy. For large distances the energy returns 
to a single-valued function of the length with approximately linear behavior. The self- 
intersection point in the energy curve indicates the presence of a first-order phase transition 
with order parameter tq. At this point we should note that exactly at the critical value, 
K^r — 63.8, the first order phase transition degenerates to a second order one. Then using 
purely thermodynamical arguments and the fact that not only the first but also the second 
derivative of the length vanishes at the critical value Tq"", we have: 

L - L„ ~ (ro - rS'^)^ (6.15) 

which is in agreement with gravity calculations using ( 16. Sp . The critical behavior vanishes 
for K, < Kcr- When the corresponding critical exponent takes the classical value 3. 

Using (16.71) and (16.151) we can calculate the behavior of the energy close to the critical point: 

Egg-E^l^ -^—{L-L,,){1-A\L-L,X^^) , (6.16) 

where A is a constant. The above critical behavior is similar to that found in [39 | B0 | [T9l [23] . 
This expression is calculated numerically by tuning k, very finely such that the first order 
phase transition, seen in L{ro), just disappears. In order to do this we must zoom into the 
region in the L{rQ) plot where both L'{rQ) and L"{rQ) simultaneously vanish. Using this 
value of K we are then able to fit the parametric plot of E and L to very high accuracy and 
obtain the classical values to high numerical precision. 

We close this section with the remark that the confining behavior appearing in the IR is 
not something expected from a gauge theory point of view. For this reason, it would be 
interesting to examine the stability of the string trajectory used for calculating the quark- 
antiquark potential under small fluctuations, through the tools obtained in [HI |12] . 

7 Discussion 

In this paper we have studied the addition of flavor to the supergravity dual of three- 
dimensional gauge theories with eight super symmetries. The unflavored background is con- 
structed by wrapping D4-branes in a two cycle of a Calabi-Yau twofold. The corresponding 
flavor branes are also D4-branes that extend along the non-compact directions of the normal 
bundle of the cycle. We have shown that the addition of these flavor branes does not break 
any of the eight supersymmetries preserved by the unflavored background. We have studied 
this system both in the quenched and unquenched approaches and we have obtained the 
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meson spectrum by analyzing the normalizable fluctuations of the probe flavor D4-branes. 
In the brane probe approximation we have also studied the Higgs branch of the theory, which 
is realized as a configuration in which the color and flavor D4-branes are recombined. 

We have succeeded in computing the backreaction of a large number of flavor branes by 
considering a continuous distribution of D4-branes smeared over their transverse angular 
directions. The flavor branes provide a source term for the RR fields that induces a violation 
of the Bianchi identity for F^. The ansatz for the backreacted background can be obtained 
by modifying the unflavored one to satisfy the new Bianchi identities. By requiring that 
the flavored background preserves the same amount of supersymmetry as the unflavored 
solution, one arrives at a system of first-order BPS equations. We have checked that the 
BPS equations imply the Maxwell-Einstein equations with sources. The functions in the 
ansatz depend on the two radial variables a and r and the BPS equations are a system 
of PDEs that must be integrated numerically. However, the solution for a = can be 
found analytically and we have shown that it reproduces the running of the coupling of the 
Yang-Mills theory with flavors. 

The supergravity solutions we have studied do not capture the rich non-perturbative struc- 
ture of the corresponding field theory duals. This fact is related to the bad IR behavior of 
the solutions, which have a naked singularity in this region. This problem is similar to the 
one presented by other similar backgrounds with the same amount of supersymmetry in 
four [HI IS] and two dimensions [22] . As in these cases, one can argue that the singularity 
can be consistently screened by an enhangon, which is the locus where the sources of the 
background become effectively tensionless and the geometry ends. The fact that the gravity 
solution reproduces the exact perturbative behavior and misses the non-perturbative effects 
seems to be due to the suppression of the latter in the 't Hooft large limit. 

We have also studied the behavior of Wilson loops in the unflavored background, which 
critically depends on the value of the parameter k of the solution. While for /t < 1/16 there 
is a maximal separation of the quark-antiquark pair, for k > 1/16 there is a transition from 
a Coulomb-like behavior at small separation to a linear potential for large separation. Above 
a certain critical value of k = k^t the transition is discontinuous (first order) while exactly at 
the critical value of n it becomes second order. At this point we have estimated the critical 
exponents, which are given by the classical mean field theory values. Actually, this critical 
behavior seems to be universal, at least for these type of models: we have found it in the 
4d system of [8|, as well as in the 2d system of [22j. In both cases there is a constant which 
is the analogue of k and a critical value for this constant. Similar results have been found 
in [391 EHl Sni [El [201 [23] in other backgrounds. Interestingly, the critical exponents in these 
systems are also given by the classical values. To properly interpret these results one should 
have a clean understanding of the meaning of the constant n (and of its analogues in other 
models). The most natural interpretation is that k parametrizes the VEV of some operator 
that is switched on. However, more work is needed to confirm this interpretation and to find 
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the precise nature of the phase transition uncovered. We are working in this direction. 
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A BPS equations 



Let us derive the BPS equations (13.121) of the general flavored background by imposing the 
preservation of eight supersymmetries. For the type of background we are studying, the 
supersymmetry variations for the dilatino A and gravitino ipM in the type IIA theory in the 
string frame are: 



5\ = 



1 

2 



^ p(4) T^MNPQ 

^ . ^1 ^ ^MNPQ J- 



Vm - 



(A.l) 



where e is a ten-dimensional spinor. The Killing spinors of the background are those for 
which 6X = 6tpM = 0. They will be characterized by a set of algebraic projection conditions 
that can be expressed in terms of products of constant Dirac matrices with fiat indices. In 
order to specify these conditions, let us choose the following vielbein basis: 



p0,l,2 



e^Z^^Rde , e^ = e*Z^i?sin 



e Z ■^do , = e *Z (dtp + cos 



(A.2) 



3^ = e *cir , = e '^r d9 



e r sm 



It is useful to express the four-form fl3.1ip in terms of the vielbein basis (lA.2p . One has: 



^4$ 



Z (^g' + \fZg j A A A + 



+ 



A e'' A A 



(A.3) 
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Then, the projections satisfied by the Kilhng spinors are: 



01234 



roi256 Til e 



(A.4) 



where the indices of the F's refer to the frame basis (]A.2|) . One can easily show that, after 
putting all these ingredients into the supersymmetry transformations for the gravitino and 
dilatino (lA.ip . we end up with the system (13.121) of first-order BPS equations for the functions 
of our ansatz. Moreover, the expression for the Killing spinors is: 



(A.5) 



where 77 is a constant spinor that satisfies the same projections (1A.4I) as e. Notice that these 
algebraic conditions imply that the system is 1/4-supersymmetric, i.e. that eight supersym- 
metries are preserved. 

Let us now show that the first-order BPS system fl3.12p implies the second-order equations 
of motion of the different fields. First of all, let us start by checking the equation of motion 
of F4, which is given by: 

d{*Fi) = 0, (A.6) 
One can prove that (1A.6P is equivalent to the following PDE: 



dr 



Zg'e' 



Z^ge 



8* -, 



a 



(A.7) 



Let us verify that equation (1A.7I) is satisfied as a consequence of the system fl3.12p . We 
will check this fact by an explicit calculation. First, by using (I3.12p . we rewrite all terms 
appearing in (]A.7p as: 



Zg'e' 



Z'^ge"" 



r2- 



cr 



ZZ 
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Z^ r 
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©(r-rg) 



R^dr 



(A.8) 



By plugging this result into ( 1A.7P one can straightforwardly verify ( ]A.6p . Thus, Fq = *F^ 
should be represented as the derivative of a five-form, dC^. Indeed, by using again (13.120 
and (lA.SP one can readily obtain an expression for this five-form potential, namely: 



= dx^ A dx^ A dx^ A 



r>2 4<I> ~ 
R Ze UJ2 



a 
Z 



da A [dtp + cos 6 dcj)) 



(A.9) 



In order to check that the Einstein and dilaton equations of motion are satisfied in the 
flavored background, after using the BPS equations, we write them in the Einstein frame. In 
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this frame we will use a one-form basis which is just the one in (1A.2P conveniently rescaled 
with the exponential of the dilaton as: 



E 



M 



e 4 e 



M 



Notice that the smearing form Q of f l3.5p can be written in this basis as: 

^ = ^^'2^''^ e'ae'ae'ae'ae'. 



(A.IO) 



(A.ii: 



Moreover, the DBI action for the smeared flavor branes in the Einstein frame takes the form: 



DBI 



-T4 f d^^x e*/V-det G^\n 

J Mio 



(A.12) 



where = e G is the Einstein frame metric and the modulus \^\^ is computed with 
G^ (from now on we will suppress the index E of G^). It follows immediately from flA.llI) 



\E . 



that \fl\ is given by: 



- 167r2 r^\Z\ ' ^^'^^^ 

The total action of the system is just the one of type IIA supergravity plus the DBI+WZ 
action of the flavor branes. The corresponding equation of motion for the dilaton is: 



1 



Again, one can show that (lA.14p is satisfied as a consequence of the system (13 .120 . The 
verification of this fact is straightforward (but tedious) and will not be detailed here. Let 
us just mention the fact that one has to evaluate the left-hand side of flA.14p by computing 
the derivatives of the first-order BPS equations in fl3.12p . In this process one generates some 
terms in which the Heaviside function 0(r — tq) is differentiated and, therefore, the Dirac 
delta function 6{r — tq) is produced. These terms match precisely the one containing Vt 
l4l) . while the remaining ones correspond to the F| term. 
The Einstein equations in the Einstein frame are: 



m 



Rmn — 2 GmnR — 2 



dM^dN^-]^GMN{d^f 
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+ T, 



MN 



(A.15) 



4(-F4 )mn - 

where Tmn is the energy-momentum tensor of the smeared flavor brane, defined as: 
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Taking into account the form of Sdbi in (1A.12P we arrive at the following expression for 
Tmn in flat components: 



Vmn " ~ 77 "1 I 
4! ^ 



(A.17) 



From this expression one readily gets the explicit values of the different components of T^tv; 
which are: 
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99 



R Nj Sir- rg) 
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T33 — T44 — T55 — Tgg 



T77 = . 



(A.18) 



By using these values one can verify, following the same strategy used to prove ( ]A.14I) . 
that the Einstein equations flA.lSP are satisfied as a consequence of the first-order equations 



B Super symmetric embeddings 

In this appendix we will characterize a class of supersymmetric embeddings of D4-branes in 
the backgrounds of sections [2] and [31 Our main tool will be kappa symmetry [13], which states 
that the supersymmetric embeddings of the D4-brane are those that satisfy the condition 
e = e, where e is a Killing spinor of the background and is a matrix which depends 
on the embedding. To write the precise form of F^, let us define the induced Dirac matrices 
on the D4-brane worldvolume as 7a = daX^^ E^^ Tj^, where X'^^C,"') are the functions that 
parameterize the embedding and are the vielbein coefficients of ten-dimensional metric. 
Then, when the worldvolume gauge field F is zero, the matrix F^ for the D4-brane is 



F. = ^ . ^ . Fn e"^-"'^ la,-as , (B.l) 
5! V-detGs 

where 'jai-as denotes the antisymmetrized product of the induced matrices, 65 is the induced 
metric on the D4-brane worldvolume and Fn is the chiral matrix in ten dimensions. Let 
us choose, as in (13.11) . x^, x^, a and ip as worldvolume coordinates and let us consider 
embeddings as in (14.111) . in which 6 and (p depend on {cr,ip) and the remaining coordinates 
are constant. Then, the kappa symmetry matrix ( IB.ll) takes the following form: 

Tk = ; = Til 72.0a,ij.2(^^ , (B.2) 

V - det G5 
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where the induced Gamma matrices are: 



7a,o,i,2 = e i 0,1,: 
7^ = i?e*Zi/2 



(9,^r3 + sin^(9,0r4 
d^6 Fs + sin 9 d^cj) T4 
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+ 



a 



Fs + cr cos6'(9o-0F6 



1 + cos 6^ 



(B.3) 



To find those embeddings that are kappa symmetric and preserve the same amount of super- 
symmetry as the original background, we should compute the action of the antisymmetrized 
product x^o-V' on the spinor. Taking into account the projections (1A.4I1 . we have: 



e ^*Fii 'jxOx^x^ai, e = [c/ + C39F39 + C49F49] e 



(B.4) 



where the coefficients appearing on the right-hand side of the above equation are: 
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d^9 da(j) - da9 d^(f) 
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da9 + cos 9 i df(f) da9 — dcr4> d^d 



(B.5) 



The right-hand side of (]B.4p should contain only the term with the identity matrix if we 
want to satisfy the F^ e = e condition for any Killing spinor of the background. Thus, we 
must demand that: 

C39 = C49 = , (B.6) 
which leads us to the following system of PDE's: 



(9w,6' — cr sin 61 da4> = 



sm 
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d(j9 + cos 9 ( d^(j) d(j9 — da-(f) 0^9 



. 



(B.7) 



The general solution of the system (IB. 71) was obtained in ref. [22]. As mentioned in subsection 
14.21 if one defines the two complex variables Ci and as in (14.121) , any holomorphic function 
Ci = /(C2) solves (IB. 71) . The corresponding embedding preserves the eight supersymmetries 
of the background. 



C Additional supergravity backgrounds 

In this appendix we find additional supergravity solutions for our unfiavored setup. First of 
all we will show that it is possible to find solutions for the unfiavored BPS system of equations 
(12. 8p that are simpler than those coming from gauged supergravity. Later on, in subsection 
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IC.ll after performing a series of duality transformations to the unflavored solution of section 
[21 we will generate a background dual to a non-relativistic system in 1+1 dimensions. 

Let us try to solve the unflavored BPS system ( 12. Sp by means of the method of separation 
of variables. Accordingly, let us adopt the following ansatz for the three functions g, $ and 
Z appearing in the background: 



9 = 9[(^) 



Z = Z^{r)Z^{a) . 



(c.i: 



It is straightforward to show that the general solution of the system (12. Sp that has the form 
of the ansatz ( IC.ll) is: 

g{a) = ci a/c4 + C3 , 
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(C.2) 



where the q's are separation constants. The ten-dimensional metric corresponding to this 
solution can be obtained by plugging the values of e~^* and Z given in (1C.2I) into a general 
ansatz. Proceeding in this way we arrive at the split 6+4 metric: 



where the six-dimensional metric is independent of cr: 

dsl = e^^dxl, + [dr^ + dnj] 
while the four-dimensional metric is independent of r: 

RV^3 
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(C.5) 

In order to rewrite the metric fIC.Sp in a more familiar form, let us perform a change of 
variables and define a new radial variable (, related to a in the following way: 



4R 



Actually, if we define a new constant a as: 



(C.6) 



(C.7) 
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the relation that gives a in terms of C, is: 



(Ci 



Clearly, C ^ which corresponds to the range a > 0. After this change of variable, one can 
easily prove that the metric fIC.SP becomes: 
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(C.9) 



which is the metric of an Eguchi-Hanson space EH^^ with resolution parameter a. When 
a = the metric (]C.9I) becomes the one corresponding to the C^/Z2 orbifold. Moreover, if 
we define the new constants t] and Q as: 



?7 = i?v^C2 
then, the warp factor e~^* becomes: 
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Using this result the six-dimensional part of the metric takes the form: 
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Notice that, for rj ^ 0, the metric flC.12p is asymptotically flat. 



C.l Non-relativistic backgrounds 



Let us now follow the procedure of ref. [45] to obtain backgrounds dual to non-relativistic 
systems by performing a combination of two T-dualities and a shift to a supergravity dual 
of a relativistic theory. First of all, we introduce light-cone variables in the standard way, 
= x° ± x^, and rewrite the initial metric and the dilaton of the unfiavored type IIA 
background as: 
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Notice that we have introduced the warp factor H. Similarly, the RR potentials in these 
variables take the form: 
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In the first step we perform a T-duality along the fiber direction ■0. In general, under T- 
duality the fields in the NSNS sector form a closed set and transform amongst themselves. 

Hence for these fields we may use the standard rules. In contrast, the transformation rules 
for the RR sector fields involve those in the NSNS sector. The metric and dilaton after this 
T-duality become: 
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while the non-vanishing NSNS and RR potentials are: 
B2 — cos 9d'^ Ad^ , 
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In the second step we perform a coordinate shift along the light-cone coordinate x~ of the 
form: 

x" ^ x~ -\- 7'0 — > dx" dx" -\- ^ydip . (C.17) 
The corresponding metric and the dilaton become: 
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while the non-vanishing NSNS and RR potentials are: 

B2 = cos 9 dtp A dcj) , 

C2 = —g^2 , = dx~^ A {dx~ + 'fdip) A dx^ A da , 
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Finally we perform another T-duality along ijj and return to a type IIA background. Now 
the metric and the dilaton become: 
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The corresponding field strengths of these potentials are Hj, = dB2, -F4 = dC^. By computing 
the exterior derivatives and using the BPS equations we have: 
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As a consistency check we can verify that these forms satisfy their equations of motion: 
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D Entanglement entropy 



In quantum field theory the entanglement entropy between two complementary spatial re- 
gions A and B is defined as the entropy seen by an observer in A who does not have access 
to the degrees of freedom of B. The holographic computation consists of finding the eight- 
dimensional surface S with minimal area such that its boundary coincides with the boundary 
of A. Then, the entanglement entropy between A and its complementary region B is given 
by the integral [16]: 

S = j^ [ d'^e-'^^s, (D.l) 

where Gio is the ten-dimensional Newton constant, given by Gio = Svr^a'^^fg and Gg is the 
induced metric on E (see also |lZj). We will consider a constant time surface S, obtained 
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by minimizing S over all surfaces that approach the boundary of A at the boundary of the 
ten-dimensional bulk manifold and that are extended along the remaining spatial directions. 
The surface we consider is not going to touch the IR region of the space, so it is feasible to 
use the UV expression ( ]2.24p . We parametrize the eight-dimensional surface in the following 
way: 

e={x,X2jA,a,e,<P,i;) , (D.2) 

with u = u{x). By computing the induced metric we end up with the following expression 
for S: 

2 7r3i?5 



S 



3G 



10 



dxu \ u'"^ + 



2* Vi?/ 



:d.3) 



Since the function S does not depend explicitly on x, the Euler-Lagrange equation derived 
from 5* can be integrated with the result being the following: 
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12 



2^2 



1/4 



:d.4) 



where Uq is the minimal value of u. From this expression we can obtain u' as a function of 
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1/2 



U 



3/2 / ^ 
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:d.5) 



Now we can compute the length / as a function of the turning point of the holographic 
coordinate: 



no 



v! 



1 



:d.6) 



We can use (ID .50 to eliminate u' in the entropy functional S. The resulting integral is 
divergent if the upper limit is infinity. For this reason we regulate it by integrating up to 
some value u^a of u and we have: 
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2 7r^R' 
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"0 
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(D.7) 



The integral appearing on the right hand side of the above equation takes the value: 



^r/2 
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(D.8) 



where B is the symbol for the incomplete Beta function. Taking the limit of the above 
expression when ^ oo we obtain finite and divergent terms: 

2 / \ 2 



HQ 



^7/2 



d^ 
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:d.9) 



Plugging the finite part of tliis result in (ID.Tp and using (1D.6P to express Uq in terms of / we 
have: 



^4 ' (D-10) 



where Rio = QsOi'^^'^ and Voo is the volume of a sphere of radius RZl^"^ along which the 
D4-branes are wrapped: 

V^=A-kR^Z,. (D.ll) 

The expression ( ID. 101) just obtained coincides with the one found in [17] for an M5-brane 
compactified in a two sphere of volume V^. This result was, of course, to be expected, given 
the UV behavior of our metric. 
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